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Abstract. In this paper, we study an version of the semiclassical approximation of magnetic 
Schrodinger operators with invariant Morse type potentials on covering spaces of compact man- 
ifolds. In particular, we are able to establish the existence of an arbitrary large number of gaps 
in the spectrum of these operators, in the semiclassical limit as the coupling constant /x goes to 
zero. 



Introduction 



A charged particle constrained to a manifold, moves along geodesic orbits. However, in the 
presence of a magnetic field, the orbit of the charged particle is no longer a geodesic, but rather 
a magnetic geodesic, that is, a trajectory determined by solving the Hamiltonian equations given 
by the symbol of the magnetic Hamiltonian. In the case of Euclidean space and constant mag- 
netic field, the magnetic geodesies are circles, and on the hyperbolic space they are the lines of 
tangency of the hyperboloids embedded in the Minkowski space. The corresponding magnetic 
Schrodinger operators are sometimes also called magnetic Hamiltonians or Landau Hamiltonians. 
The magnetic Schrodinger operators turn out to possess rich spectral properties cf. 



, 1 16], that are important in the analysis of the quantum Hall effect cf. H 




We begin by reviewing the construction of the magnetic Hamiltonian on the universal covering 
space M of a compact connected manifold M. The fundamental group F = 7ri(M) acts on M 
by the deck transformations, so that M/F = M. Let u; be a closed 2-form on M and B be its 
lift to M, so that i? is a F-invariant closed 2-form on M. We assume that B is exact. Pick a 
1-form ^ on M such that dA = B. As in geometric quantization we may regard A as defining 
a Hermitian connection V = d + iA on the trivial line bundle C over M, whose curvature is iB. 
Physically we can think of A as the electromagnetic vector potential for a magnetic field B. Using 
the Riemannian metric the magnetic Laplacian is given by 

Ha = V*V = {d + iAy{d + iA). 
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Then the Schrodinger equation describing the quantum mechanics of a single electron of mass m 
which is confined to move on the covering space M in the presence of a periodic magnetic field is 
given by 

d 1 

ih—il^ = {hd + ieA)* {M + ieA) i) + ii~'^Vip 

where Ha is the magnetic Laplacian, V \s a. F-invariant electric potential function, h is Planck's 
constant, e is the electric charge of the charge carrier and /i is the coupling constant. In the time 
independent framework, the relevant operator is the magnetic Schrodinger operator 

where the physical constants are set equal to 1. It is the qualitative aspects of the spectrum 
of HAyifj) that are relevant to the study of the quantum Hall effect. An important feature of 
the magnetic Schrodinger operator is that it commutes with magnetic translations, that is with a 
projective action of the fundamental group F. Let a denote the multiplier or U (l)-valued 2-cocycle 
on F defining this projective action. Under the assumption that the Kadison constant of C*(F, a) 
is positive, it was proved by Briining and Sunada 0] that there are only a finite number of gaps in 
the spectrum of the magnetic Schrodinger operator HAy that lie in any left half-line (-co. A]. In 
fact, under the same hypotheses, they obtain Weyl-type asymptotics of the number of gaps in the 
spectrum, as A — > oo. On the other hand, it was proved in [^] that the Kadison constant is positive 
under the assumptions that the Baum-Connes conjecture with coefficients holds for F, and that 
F has finite cohomological dimension and finally that a defines a rational cohomology class. The 
assumptions on F are satisfied for instance when F is a discrete subgroup of SO{l,n) or SU{l,n) 
or of an amenable Lie group. Our paper on the other hand proves that HAy can have arbitrarily 
large number of gaps, if V is chosen to be a suitable Morse-type potential. An outstanding open 
problem is to construct magnetic Schrodinger operators HAy that have infinitely many gaps in 
some half-line (— oo. A], for any F, where of course a has to be an irrational cohomology class. This 
has been shown to be the case at least in the discretization of the magnetic Schrodinger operator, 
called the Harper operator, and when F = Z^, cf. 0,0, There have been other interesting 
papers on gaps in the spectrum of magnetic Schrodinger operators on covering spaces. One such is 
1 14 1 where the asymptotics of the size of gaps in the spectrum are obtained on hyperbolic spaces. 



as the mid-point of the gap tends to infinity, generalizing known results on Euclidean spaces. In 
two dimensions, magnetic Schrodinger operators on covering spaces are the Hamiltonians in the 
model for the integer quantum Hall effect when the covering space is the Euclidean plane S, and 
for the fractional quantum Hall effect when the covering space is the hyperbolic plane [|^, Q. 
The standard assumption made is that the Fermi energy level lies in a gap of the spectrum of the 
Hamiltonian (which however can be relaxed to the assumption that it lies in a gap in extended 
states of the Hamiltonian after further analysis). 

In this paper, we study the semiclassical approximation, as the coupling constant approaches 
zero, of the Hamiltonians HAyifJ-), for a class of Morse type potentials V, which include all 

functions V = \df\'^, where / is a F-invariant Morse function on M, that is / is the lift to M of a 
Morse function on M. We show that the spectrum of this operator, is approximated by the union 
of the spectra of model operators which are defined near the critical points of the Morse potential 
V. In particular, we are able to deduce the existence of an arbitrarily large number of gaps in the 
spectrum of these Hamiltonians, for n sufficiently small, a fact which is of crucial importance in 
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the study of the quantum Hall effect. We adapt the Lp' version of semiclassical approximation of 
Witten m, @, [||] to our context. 

The paper is organized as follows. We first give a summary of our main results. Then we recall 
some preliminary material on projective unitary representations and the associated von Neumann 
algebra of operators commuting with this algebra, together with the von Neumann trace and 
von Neumann dimension function. In the next section, we establish the first main theorem on 
the existence of spectral gaps for large values of the coupling parameter or equivalently for small 
values of the coupling constant /i. Here we use the method of semiclassical approximation. 



Summary of main results. In this paper, we will assume that the potential y is a smooth 
T-invariant function which satisfies the following Morse type condition: V{x) > for all x £ M. 
Also if V{xo) = for some xq in M, then there is a positive constant c such that V{x) > c|x — xqP 
for all X in a neighborhood of xq. We will also assume that V has at least one zero point. 

We remark that all functions V = where \df\ denotes the pointwise norm of the differential 
of a F-invariant Morse function / on M, are examples of Morse type potentials. 

We next enunciate the principle on which this paper is based. Associated to each Hamiltonian 
Ha,v{i^)i there is a model operator K (cf. section 2) which has a discrete spectrum. It is defined 
as a direct sum of harmonic oscillators, associated with the potential wells of y in a fundamental 
domain of F in M. 

Semiclassical approximation principle: Let V he a Morse type potential. Then in the semi- 
classical limit as the coupling constant ^ goes to zero, the spectrum of the Hamiltonian HAyipi) 
"tears up" into hands which are located near the eigenvalues of the associated model operator K. 

The following main theorem establishes the existence of arbitrarily large number of gaps in the 
spectrum of HAyifJ-) whenever ^ is sufficiently small. The proof of this theorem uses an analogue 
of Witten's semiclassical approximation technique for proving the Morse inequalities. The L^- 
analogue (in the absence of a magnetic field) was proved by Shubin [^] , also Burghelea et al. |5| . 



We modify the proof in |19] to obtain the result in section 2. The physical explanation for the 
appearance of gaps in the spectrum HAyifJ-) is that the potential wells get deeper as — > and 
the atoms get (asymptotically) isolated, so that the energy levels of Ha,v{ij) are approximated 
by those of the corresponding model operator K. 

Theorem (Existence of spectral gaps). Let V he a Morse type potential. If E € C is such that 
E ^ spec(K), then there exists no > such that for all ^ G E is in the resolvent set of 

HAyip)- If in addition E is real and lies between two eigenvalues of K, then E is in a spectral 
gap of HAyifJ-). 

Since the spacing between the eigenvalues of K is bounded below, it follows that there exists 
arbitrarily large number of gaps in the spectrum of Ha.v{^j) whenever the coupling constant fi is 
sufficiently small. 

Remark. We observe that the operator Ha + fJ-~^V = ^~^Ha.v{p) also has arbitrarily large 
number of gaps in its spectrum whenever the coupling constant /i is sufficiently small. Analogous 
results in the special case of Euclidean space were obtained in [p^] . 
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1. Preliminaries 



Let M be a compact connected Riemannian manifold, F be its fundamental group and M be its 
universal cover, i.e. one has the principal bundle T ^ M ^ M. To make the paper self-contained, 
we include preliminary material, some of which may not be new, cf. [HI, P], |15], [16|. 



1.1. Projective action, or magnetic translations. Let w be a closed real- valued 2-form on 
M such that B = p*uj is exact. So B = dA where A is a 1-form on M. We will assume A without 
loss of generality that A is real- valued too. Define V = d+ iA. Then V is a Hermitian connection 
on the trivial line bundle over M with the curvature (V)^ = iB. The connection V defines a 
projective action of T on functions as follows. 

Observe that since B is P-invariant, one has = j*B — B = d{j*A — A) V7 € P. So 7*^ — ^ 
is a closed 1-form on the simply connected manifold M, therefore 

YA - A = d^j^, V7GP, 

where is a smooth function on M. It is defined up to an additive constant, so we can assume 
in addition that it satisfies the following normalization condition: 

• tp^{xo) = for a fixed xq £ M, V7 € P. 

It follows that tp^ is real- valued and ipeix) = 0, where e denotes the neutral element of P. It is 
also easy to check that 

• ip-yix) + ipyi'yx) — ipY^{x) is independent of x G M, V7,7' € P. 

Then (7(7,7') = 6xp(— Z'0^(7' • xq)) defines a multiplier on P i.e. a : P x P ^ U{1) satisfies 

• cj(7,e) = a(e,7) = 1, V 7 G P; 

• o-(7i,72)o-(7i72,73) = o-(7i,7273)c^(72,73), V7i,72,73 G P (the cocycle relation). 
It follows from these relations that (7(7,7""'^) = (7(7"^, 7). 

The complex conjugate multiplier (1(7,7') = exp(i'0^(7' • xq)) also satisfies the same relations. 

For u € L'^{M) and 7 € F define 

U^/U = (7^"^)*n, S^u = exp(— i'i/'7) u. 
Then the operators = U^ o S-y satisfy 

Te = Id, T^iT^a = f^(7l,72)7'7i72; 

for all 71,72 G F- In this case one says that the map T : F — > U{L'^{M)), 7 1-^ T^, is a projective 
(P, (j)-unitary representation, where for any Hilbert space Ti we denote by hl{T-L) the group of all 
unitary operators in 7i. In other words one says that the map 7 1— > defines a {T,a)-action 
in n. 
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It is also easy to check that the adjoint operator to in L^(M) (with respect to a smooth 
F-invariant measure) is 

The operators are also called magnetic translations. 

1.2. Twisted group algebras. Denote by ^^(F) the standard Hilbert space of complex-valued 
L^-functions on the discrete group F. We will use a left (F, fT)-action on ^^(F) (or, equivalently, a 
(F, a)-unitary representation in ^^(F)) which is given explicitly by 

Ti;f{i) = /(7~S')^(7, 7-S'), 7, 7 e r. 

It is easy to see that this is indeed a (F, fT)-action, i.e. 

= Id and T^^^T^^ = a{^i , ^2)T^,^, , V71 , 72 e T. 

Also 

(T^)* = a(7,7-')rili. 



Let 



A^{T, a) = {Ae B{e\r)) : [T^, A] =0, V7 e f} 



be the commutant of the left (F, a)-action on £^(F). Here by 13{Ti.) we denote the algebra of all 
bounded linear operators in a Hilbert space Ti. By the general theory, ^^(F, a) is a von Neumann 
algebra and is known as the (right) twisted group von Neumann algebra. It can also be realized 
as follows. Let us define the following operators in ^■^(F): 

T^''/(7') = /(7'7y(7',7), 7,7' er. 

It is easy to check that they form a right (F, cj)-action in ^^(F) i.e. 

r,« = Id and T^J^^ = a{-fi,j2)T^^^^, V7i,72Gr, 

and also 

(r^«)* = a(7,7"')T,^-i. 
This action commutes with the left (F, a')-action defined above i.e. 

T^Tf = T:^T^, V7,7'gF. 

It can be shown that the von Neumann algebra A^{T,a) is generated by the operators {T^j^gr 
(see e.g. a similar argument in flF 



Similarly we can introduce a von Neumann algebra 

A^{r,a) = [Ae B{e\r)) [t^,a] = o, V7 g f}. 

We will refer to it as (left) twisted group von Neumann algebra. It is generated by the operators 
{T^l-ygrj and it is the commutant of A^{T,a). 
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Let US define a twisted group algebra C(r,cr) which consists of complex valued functions with 
finite support on T and with the twisted convolution operation 

{f*9){l)= XI /(7i)S'(72)c^(7i,72)- 

71:72:7172 =7 

The basis of C{T,a) as a vector space is formed by 5-functions {d-yj-^er, ^-yil') = 1 if 7 = 7' and 
otherwise. We have 

'^71 * '^72 = f^(7l,72)'57i72- 

Associativity of this multiplication is equivalent to the cocycle condition. 

Note also that the J-functions {(5-y}^gr form an orthonormal basis in ^^(F). It is easy to check 
that 

T^Sy = ct(7, 7')V' ^7^'^7' = <t(7'7"\ 7)'^7'7-i- 

It is clear that the correspondences 6^ 1— > and (5-y 1— > define representations of C(r,fT) 
and C(r, a) respectively. In both cases the weak closure of the image of the twisted group algebra 
coincides with the corresponding von Neumann algebra {A^{r, a) and A^{T, a) respectively). The 
corresponding norm closures are so called reduced twisted group C*- algebras which are denoted 
C*{T,a) and C*{r,a) respectively. 

The von Neumann algebras A^{T,a) and A^{r,a) can be described in terms of the matrix 
elements. For any A E B{£'^{T)) denote Aa^p = {A6f3,6a) (which is a matrix element of A). Then 
repeating standard arguments (given in a similar situation e.g. in |l8|) we can prove that for any 
A G B{£'^{r)) the inclusion A G A^{T,a) is equivalent to the relations 

Ayxay = ^(7, a;)cr(7, y)A^,.y , Vx, y, 7 G T. 

In particular, we have for any A G A^(r, a) 

■'^"/x,'yx — ■'^x,x 5 Vx,7 £ r. 

Similarly, for any A G B{£'^(r)) the inclusion A G A^{T,a) is equivalent to the relations 

— a{x,^)a{y,j)Ax^y , Vx,y, 7Gr. 

In particular, we have 

for any A G A^{T,a). 

A finite von Neumann trace trr^^ '■ ^^(r,^) — > C is defined by the formula 

tTT,aA= {A6e,6e). 

We can also write trr,^ A = A^^^ = {A5-y, (5^) for any 7 G F because the right hand side does not 
depend of 7. 

A finite von Neumann trace trr^o- : A^{T,a) ^ C is defined by the same formula, so we will 
denote by trp any of these traces. 

Let Ti. denote an infinite dimensional complex Hilbert space. Then the Hilbert tensor product 
£^{T) is both (F, CT)-module and (F, cT)-module under the actions 7 1-^ Ti" 1 and 7 1— > TI^ (S) 1 
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respectively. Let Ay^{T, a) and yl:^(r, a) denote the von Neumann algebras in £'^{T)(E)'H which are 
commutants of the {T,a)- and (F, CT)-actions respectively. Clearly A^{T,a) = A^{T,a) B{H) 
and ^^(r, a) ^ A^{T, a) (g) B{n) in the usual sense of von Neumann algebra tensor products. 

Define the semifinite tensor product trace Trr = trp CS) Tr on each of the algebras A^{T, a) and 
A^{T,a). Here Tr denotes the standard (semi-finite) trace on B{H). 

Let us recall that a closed linear subspace V C TC is affiliated to a von Neumann algebra A of 
operators in a Hilbert space TC if Py S A where Py is the orthogonal projection on V in 7i. This 
is equivalent to saying that V is invariant under the commutant A' of A in 7Y. 

The von Neumann dimension dimp of a closed subspace V of ^^(F) 0Ti. that is invariant under 
{T^ (8) 1, 7 € r} (or, equivalently, affiliated to A^{T, a)) is defined as 

dimr(F) = Trr(Py). 

The same formula is used for the subspaces which are invariant under {TJ^ (g) 1, 7 G F} (or, 
equivalently, affiliated to Ay^{T,a)). 

Also the von Neumann rank of an operator Q G Ay^{T,a) or , cj) is defined as 

rankr((5) = dimr(Range((5)), 
where the bar over Range((5) means closure. 

1.3. Magnetic Hamiltonians. The magnetic Laplacian on L?'{M) is defined as 

Ha = V*V = {d + iA)* {d + iA) 
and more generally, the magnetic Hamiltonian or Magnetic Schrddinger operator is defined as 

where V is any F-invariant smooth function on M. The Hamiltonian H = HAyifj) is a self 
adjoint second order elliptic differential operator. It commutes with the magnetic translations 
(for all 7 G F), i.e. with the (F, cj)-action which was defined above. To see this note first that the 
operators C/^ = (7^^)* and (the multiplication by exp(— z-i/;^)) are defined not only on scalar 
functions but also on 1-forms (and actually on p-forms for any p > 0) on M. Hence the magnetic 
translations are well defined on forms as well. The operators are obviously unitary on the 
spaces of forms, where the LP' structure is defined by the fixed F-invariant metric on M. An easy 
calculation shows that T^V = VTy on scalar functions. By taking adjoint operators we obtain 
T^V* = V*Ty on 1-forms. Therefore T^Ha = HaT^ on functions. Since obviously T^V = VT^, 
we see that HAyifJ-) commutes with T-y for all 7. 

In dimension 2, it is the spectrum and the spectral projections of the magnetic Schrodinger 
operator that are of fundamental importance to the study of the quantum Hall effect. We remark 
that it is virtually impossible to explicitly compute the spectrum of Ha.v for arbitrary V which is 
F-invariant, even in 2 dimensions and for simplest manifolds. Nevertheless, Comtet and Houston 
1^] computed the spectrum of Ha (with y = 0) on the hyperbolic plane with the magnetic 
potential A = 0y~^dx (which corresponds to the constant magnetic field), where we can assume 
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without loss of generality that ^ > 0. This spectrum is the union of a finite number of eigenvalues 
{{2k + 1) 9 - k{k + 1) : k = 0,1,2, - ■ ■ < O - ^} and the continuous spectrum [i + 9^, oo) . 

Since H commutes with the (F, (T)-action, it follows by the spectral mapping theorem that the 

spectral projections of H, Ex = X(-oo,x]{H) are bounded operators on L'^{M) that also commute 
with the (r, cr)-action i.e. T^Ex = ExT^, V 7 G F. The commutant of the (F, (T)-action is a von 
Neumann algebra 

Uj^ir,a) = {q € BiL\M)) : T^Q = QT^, V 7 e f} . 

To characterize the Schwartz kernels kQ^x, y) of the operators Q G Uj^iV, a) note that the relation 
T^Q = QT^ can be rewritten in the form 

eiV'7WA;Q(7x,7y)e-^'^^(2') = kQ{x,y), 'ix,y eM V7 G F, 

so Q G hlj^j{T,(f) if and only if this holds for all 7 G F. In particular, in this case kQ(x,x) is 
F-invariant. For the spectral projections of H we also have Ex G Uj^jiV, a), so the corresponding 
Schwartz kernels also satisfy the relations above. Note that the Schwartz kernels of Ex are in 
C~(M X M). 

To define a natural trace on U^{T,a) we will construct an isomorphism of this algebra with 
the von Neumann algebra A}f^{T, a), where H = LP'{T). 

Let ^ be a fundamental domain for the F-action on M. By choosing a connected fundamental 
domain T for the action of F on M, we can define a (F, cr)-equivariant isometry 

(1.1) \5 ■.L^(M)'^f{T)(^L^{T) 

as follows. Let i : ^ M denote the inclusion map. Define 

U(0) = ^5^® i*{T^(f>), (f> e L^iM). 

Lemma 1.1. The map U : L^{M) —>■ £^(F) (g) L^{J-') defined above is a {y,a)-equivariant 
isometry, where the {T,a)-action is given by the operators T-y and (g) 1 on the spaces L'^{M) 
and i'^{r) (g) LP'{J^) respectively. 

Proof. Given (f) G L'^{M), we compute 

U(r^0) = ^ <5y (2) i*{TyTj(l)) = ^ (T(7',7)(5y (g i*{T^>^^) 
yer yer 

yer 

which proves that U is a (F, c7)-equivariant map. It is straightforward to check that that U is an 
isometry. □ 
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Since Wj^(r, u) is the commutant of {Xy}-ygr, and A^iJ^, a), with H = L^{J^), is the commutant 
of {r^ (8) l}7gr, we see that U induces an isomorphism of von Neumann algebras U^{r,a) and 
A^{T,a). Therefore we can transfer the trace Trr from Ay^{r, a) to U^{T, a). The result will be 
a semifinite T-trace on Uj^{T,a) which we will still denote Trp. 

It is easy to check that for any Q G U-^{r,a) with a finite F-trace and a continuous Schwartz 
kernel kg we have 

Trr Q = J kQ{x,x)dx 

where dx means the F-invariant measure. An important particular case is a spectral projection 
Ex of the magnetic Schrodinger operator H = Ha y as considered above. The projection E\ has 
a finite F-trace and a C°° Schwartz kernel. Therefore we can define a spectral density function 

Nr{X;H) = TrrEx, 

which is finite for all A G R. It is easy to see that A h- > Ny{X;H) is a non-decreasing function, 
and the spectrum of H can be reconstructed as the set of its points of growth, i.e. 

spec(i7) = {A G M : Arr(A + e;H)- Nr{X -e;H)> 0, Ve > 0}. 

The von Neumann dimension dimp of a closed subspace V of L?[M) is well defined if V is 
invariant under T^, V7 G F. Also the von Neumann rank of an operator Q G W^(F,a) is well 
defined. 



1.4. Gauge invariance. If we make another choice of vector potential A! such that dA! = B, 
then it follows that ^' — ^ is a closed 1-form on a simply connected manifold M, and therefore 
it is exact, i.e. A' = A + dcj), where cf) G C°°{M). We will always assume that (/) is normalized by 
the condition ^(xo) = 0. It follows that is real- valued. 

It then follows that the connection V = d + iA gets unitarily conjugated into a new connection 
V' = d + iA' = e-^'^Ve''^. Therefore V* = e-^'^V*e^'^ and H' = e'^'t' H e^'t' , where H = Ha,v, 
H' = i?^/ y,and V is F-invariant. In particular, H' and H are unitarily equivalent. Let us repeat 
the constructions of the previous subsections with A replaced by A' indicating relations of the 
modified objects with the old ones. 

7* A' — A' , il>iy{xo) = (with the same point G M as 

o"' (71,72) = exp (-^■0^^(72x0)) = (7(7i,72)exp (-i [(^(7172X0) - (l>{'yiXo) - <?!)(72a;o)]) ■ 

The modified magnetic translations are defined by = U'^S'^, where U!^ = = (7^^)* and 
S' = exp(-z^;). Then T'T' = a'{^i,j2)T' 



Define the function i/j^ from dtp^ 
above). Then 

The new cocycle will be 
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Clearly H' commutes with the modified magnetic translations T^, V7 € T. The relation between 
old and new magnetic translations is 

which is again the same unitary conjugation up to a constant unitary factor. 
Now we can introduce a von Neumann algebra 

ZY^(r,a') = {Q' G B{L\M)) : T'^Q' = Q'^, V7 G T}. 

Clearly the map 

Z^^(r, a) Uj^iT, a'), Q ^ Q' = e-'^Qe't 

is an isometric isomorphism of von Neumann algebras. It is easy to see that this isomorphism 
preserves the F-trace which is defined on both algebras. If an operator Q has a smooth Schwartz 
kernel kQ{x,y) and a finite F-trace, then 

hence 

kQ>{x,x) = kQ{x,x), Vx G M, 
so the equality Trr Q' = Trr Q follows from the expression of the F-traces in terms of kernels. 

If Ex and E'^ are spectral projections of H and H' = e~^^He^'^ respectively, then clearly 
E'^ = e-'-'t'Exe'-'t', therefore Trp^^^ = Trr^A- This means that the spectral density functions of H 
and H' coincide, i.e. the spectral density function is gauge invariant. 



2. Semiclassical approximation and the existence of spectral gaps 



We will study an L^-version of semiclassical approximation, which is similar to the ones which 
appear when we take the Witten deformation of the de Rham complex and consider the corre- 



sponding Laplacian (cf. pO| , 10, 12]. For the case of the algebra corresponding to the regular 
representation of 7ri(M), such asymptotics were first proved in [^], see also |5| for a related semi- 
classical approximation technique. The proofs given in |jl^ will be adapted to work in the more 
general situation that we need in this section. 

Recall that H = H^y^^) = + ^x^^V is a second order differential operator acting on 

L'^{M), such that it commutes with the projective unitary (F, o")-action on L'^(M), given by the 
magnetic translations {T^, 7 G F}. Moreover, note that H is a second order elliptic operator with 
a positive principal symbol, order operator, is a non-negative potential function on M which 
has only nondegenerate zeroes and ^ > is a small parameter. 

Actually the results will not change if we add to H any F-invariant zeroth order operator, i.e. 
multiplication by a smooth F-invariant function. 

Let us recall that Ny{X] H) denote the von Neumann spectral density function of the operator 
H which can be defined as 

Nt{\H) = TttEx{H) = Trr {xi-oc.,x]{H)) , 
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where xf means the characteristic function of a subset F cM. 

Let us choose a fundamental domain ^ C M so that there is no zeros of V on the boundary of 
T. This is equivalent to saying that the translations {7.F, 7 € T} cover the set y-i(O) (the set 
of all zeros of V). Let ¥"^{0) CMF = {xj\j = 1, . . . , be the set of all zeros of V in JT; Xi 7^ Xj 
if i j. Let K denote the model operator of H (cf. |jl9|), which is obtained as a direct sum of 
quadratic parts of H in all points xi, . . . , xat. More precisely, 

where Kj is an unbounded self-adjoint operator in L'^(W^) which corresponds to the zero xj. It is 
a quantum harmonic oscillator and has a discrete spectrum. We assume that we have fixed local 
coordinates on M in a small neighborhood B(xj,r) of xj for every j = 1, . . . ,N. Then Kj has 
the form 

where all the components are obtained from H as follows. In the fixed local coordinates on M 

(2) 

near xj, the second order term Hj is a homogeneous second order differential operator with 
constant coefficients (without lower order terms) obtained by isolating the second order terms 

in the operator H and freezing the coefficients of this operator at xj. (Note that Hj does not 

(2) 

depend of A.) The zeroth order term Vj is obtained by taking the quadratic part of V in the 
chosen coordinates near Xj. 

More explicitly, 

i,fc=l * i,k=l * 

where {g^^) is the inverse matrix to the matrix of the Riemannian tensor {gik)- 



We will say that H is flat near xj \i H = Kj near Xj. (In particular, in this case we should 
ve j4 = near Xj.) 

We will also need the operator 



K{ii) = ®i<j<NKj{ij), 

where 

K,{ij) = ^,Hf + H~^vf\ /x>0. 
It is easy to see that K[ii) has the same spectrum as K = K{\). 

Let {op : p G N} denote the set of all eigenvalues of the model operator K, 7^ aq for 
p ^ q, and rp denote the multiplicity of Op, i.e. rp = dimKer [K — Upl) where K is considered in 
{Lp'iW^))^ . Denote by N{X;K) the distribution function of the eigenvalues of K, i.e. N{X;K) is 
the number of eigenvalues which are < A (multiplicities counted). 

The following is the main result of this paper. 
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Theorem 2.1 (Semiclassical Approximation). For any R > there exist constants C > and 
fiQ > such that for any fi € (0, /Uq) and X < R, one has 

N{X - C/x^/^ K) < Nr{X; H) < N{X + Cfi^'^; K). 

Therefore 

oo 

spec(iJ) n (-00, -R] C y (op - Ci//^, ap + C/x^/^) , 
p=i 

where {ap : p € N} denotes the set of all eigenvalues of the model operator K. Moreover for any 
p = 1, 2, 3, . . . with Op G [—R, R] and any fi G (0, fio) one has 

Nr{ap + Cfi^/^; H) - Nr{ap - Cfi^^^; H) = rp = N{ap + 0; - N{ap -0; K). 

This means that for small values of /i, the spectrum of H concentrates near the eigenvalues of 
the model operator K, and for every such eigenvalue, the von Neumann dimension of the spectral 
subspace of the operator H, corresponding to the part of the spectra near the eigenvalue, is exactly 
equal to the usual multiplicity of this eigenvalue of K. 

Proof of Theorem on Existence of Spectral Gaps. Let us assume that A € M and A ^ spec(-ftr). 



and then choose R > X. Let C > and /xq > be as in Theorem |2.1| . Then by taking even 
smaller fiQ we will get 

oo 

p=i 

Therefore A spec (ff). □ 

The proof of Theorem 2.1 will be divided into 2 parts: estimating Ny'(X,H) from below and 
from above. 

2.1. Estimate from below. We will start by proving an estimate from below for Nr{X;H) 
where H = H{^), X < R with an arbitrarily fixed i? > 0, and /x J. 0. 

Our main tool will be the standard variational principle for the spectral density function. 
Lemma 2.2 (Variational principle). For every A G M 

(1) iVr(A;//)=sup{dimrL|LcDom(i/); {H f, f) < X{f, f), yfeL}. 

It is understood here that L should be a closed (T, a) -invariant subspace in L'^(M), i.e. closed 
subspace which is invariant under all magnetic translations Ty, 7 G L. 



A similar variational principle for the usual action of T was used e.g. in |11, 19|. 
We will now describe an appropriate construction of a test space L. 

Fix a function J G C^(M") such that < J < 1, J(x) = 1 if |x| < 1, J{x) = if \x\ > 2, and 
(1 — J^)^/^ G C°°(M"). Let us fix a number k, < k < 1/2, which we shall choose later. For any 
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fi > define J^^\x) = J(/i^'^x). This will be our standard cut-off function. Let Jj = J^^^ in the 
fixed coordinates near Xj. Denote also Jj^^ = (7^^)*Jj. (This function is supported near 7%.) 

We will always take fi G (0, /xq) where no is sufficiently small, so in particular the supports of 
all functions Jj^^ are disjoint. Denote 



1/2 

hi 



Clearly, Jq G C°°(M). Note that there exists cq > such that V > cofi^'^ on supp Jq. 

functions Jj^-y operators. 

Denote by {ipm,j \ = !> 2, . . . } an orthonormal system of eigenfunctions of the operator Kj 
in L^(R") where the coordinates in a neighborhood of the origin E M" are identified with the 
chosen coordinates near Xj, so that corresponds to Xj. The corresponding eigenvalues will be 
denoted \m,j- 

Let us define (pmj = Jji^mj, extended by outside of a fixed small ball centered at xj. Then 
(f'mj G Cq°{M) and it is supported near xj. 

Lemma 2.3. If 1/3 < k < 1/2, then the functions (j)m,j satisfy the following "almost orthogonal- 
ity" relations 

{4>m,j,4'm',j') = 5j,j'{5m,m' + 0{fl'^)), 
{H4>m,j,<t>m',j') = Sj,j'{^m,jSm,m' + 0{lJ?'^ ^)), 

where j = 1, . . . , N , and m belongs to a finite set. 



The proof is the same as the proof of Lemma 2.3 in [19|. The unboundedness of A does not 
matter because only a finite number of points xj is involved. 

Not define 4>m,j,'y = T-yfpmj- Then 4>m,j,'r is supported near jxj. 

Lemma 2.4. If 1/3 < k < 1/2, then the functions (pmj satisfy the following "almost orthogonal- 
ity" relations 

{II(l>m,j,-y,(t>m',j',Y) = ^ {^mj^m,m' + 0{fl^'^^^)), 

where j = 1, . . . , A^, 7 G F, and m belongs to a finite set. 

Proof. The first relation is obvious because the operator Ty is unitary and moves supports by the 
action of 7. To prove the second estimate note that it is obvious if 7 7^ 7'. If 7 = 7', then we get 

{II4>m,j,j,(t>m',j',j) = {IIT-y(t)m,j,Tj4'm.',j') = {T^^ HTj(f>mJ , (pm' ,j') = {H 4'm,j , 4>m' : 

and we can use the previous Lemma. □ 



14 V. MATHAI AND M. SHUBIN 

Now we will define two closed linear subspaces in C L^i^) and $a <^ L'^{M) as follows: 

= span{(pm,j\ >^m,j < A}, 

$A = spaif{(l)m,j,'y\ \m,j < A}, 

where span^ stands for closed linear span. Clearly dim$^ < oo, and $a is a closed (F, c7)-invariant 
subspace in L^(M). 

Lemma 2.5. dim${^ = dimr $a = -^(A; i^)- 

Proof. Clearly = ©^Li'^'A^j) where is spanned by {(pmj} with fixed j. But we have 
dim$^^. = N{\;Kj). Indeed, the eigenfunctions {ijjm,j\m = 1, . . . , N{X; Kj)} are linearly in- 
dependent and real analytic, so the corresponding 4>m,j are also linearly independent because 
<pm,j = i'mj near xj. It follows that dim$^ = N{X;K). 

For any fixed j denote by {(j)m,j \ "m = 1, • • • , -^'^(A; Kj)} the orthonormal system which is ob- 
tained from the system {(pmjlm = !)•••) -^(A; i^j)} by the Gram-Schmidt orthogonalization 
process. Then {(t>m,j \ ^mj < A} is an orthonormal basis in and {T^cpmjl A^j < A,7 G F} is 
an orthonormal basis in Note that T^cpmj € Cq°{M) and it is supported near ^Xj. 

Denote by and Px the orthogonal projections on and <I>a respectively, and Kx their 
Schwartz kernels. Then 

Kx{x,y)= ^ 4'm,j®4'm,j= ^ ^rn,j{x)(i)m,j{y)- 

{m,j\\m,j<>-} {m,j\\rn,3<>'} 

In particular, 

{m,3\ A^j<A} 

Similarly we find 

Kx{x,y)= E Ty4>m,j ^ T^4>m,j 

{m,j,7| Amj <A} 

and 

{m,j,7l Amj<A} 

It follows that Kx{x,x) = K^{x,x) for all x e J-', therefore 

dimr $A = Trr(PA) = J Kx{x, x)dx = J K^{x, x)dx = Tr(Pf ) = dim $f = N{X; K), 
which proves the lemma. □ 

Proposition 2.6. For any R > and k G (0, 1/2) there exist no > and C > such that for 
any X < R and any ji G (0, /Uq) 

(2) Nr{X + CiJi^-H)>N{X-K). 
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Proof. Note first that if the estimate (^) holds with some k > 0, it holds also for all smaller values 
of K. Now we should argue as in the proof of Lemma 2.6 in ||l^ to conclude that <I>;^ G Dom{H) 
and 

{Hf,f)<{X + Cfi^^-'){f,f), /G$A. 

Applying the variational principle (Lemma |2 .21 ) we conclude that the estimate (|2|) holds with 3k— 1 
instead of k. It remains to notice that 3k— 1 takes all values from (0, 1/2) when 1/3 < k < 1/2. □ 



Corollary 2.7. For any R > and k € (0, 1/2) there exist /xq > and C > such that for any 
X < R and any n G (0, fio) 

(3) Nr{X;H) > N{X-Cfi^;K). 



Remark 2.8. If H is flat near each of the points xj, j = 1, . . . ,N, then a better estimate is 
possible. Namely, for any i? > 0, k E (0, 1/2) and e > there exist > and C > such that 
for any X < R and any /x G (0, /io) 

(4) Nt{X;H) > N{X - C exp{-C-^fi-^+'y,K). 

(Arguments in Sect. 2 of |[l^ apply here too.) 

2.2. Estimate from above. Here we will prove an estimate from above for Ny{X;H), similar 
to the estimate (Q). 

Lemma 2.9. Assume that there exists an operator D £ U-^{T,a) (i.e. a bounded operator in 
LP'{M) commuting with all magnetic translations T^, 7 G Lj, such that 

(a) rankr D < k; 

(h) H + D>X. 

Then Nr{X - e; H) < k for any e > 0. 



The proof does not differ from the proof of Lemma 3.7 in | |19| . 

We will prove that for any i? > there exist C > and /xq > such that for any X < R and 
/i G (0,//o) there exists an operator D G U^{T,a), which satisfies the conditions (a), (b) with 
k = N{X; K),X = X- Cfi^/^. It would follow from Lemma |]| that 

(5) NriX-C^^/^;H) < N{X;K), 
hence 

(6) Nr{X;H) <N{X + Cii^^^;K), 
which is the desired estimate from above. 

Denote e{^\h) — X{-oo,x]{^j{f^)) which is the spectral projection of Kj^fj,). It is an operator 
of the finite rank N{X;Kj) in L^(M"). Identifying the coordinates in a neighborhood of G M"" 
with the chosen coordinates near Xj, so that G M" corresponds to xj, we can form an operator 
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Dj = LJjE^\fi)Jj in L'^{M), where L G M, L > A. Clearly Dj has a smooth Schwartz kernel 
supported in a neighborhood of {xj,Xj) in M x M. 

Now denote 

n 

and 

D = Y,T^DrT-\ 

It is easy to check that T^D = DT^^ V7 G F. 

Note that rank£^^-'^(/u) = N{X;Kj), hence rankDj < N{X;Kj) and rankZ)^ < N{X;K). (In 
fact, it is easy to see that we have equalities in both inequalities above, but we do not need this.) 
Also Djr has a Schwartz kernel in C^{J- x T). 

Denote by Pjr the orthogonal projection on the image of Djr. Choosing an orthonormal basis 
in this image we see that the Schwartz kernel of Pjr is also in Cq°{J- x J^). For the orthogonal 
projection P on the closure of the image of D we have 

Clearly rank Djr = Tr Pjr and rankp D = Trr P. But the argument from the proof of Lemma p.5| 
of the previous subsection shows that Trr P = Tr^jr, so 

rankrl) = rankZ)^ < N{X;K), 

which proves the condition (a) with k = N(X; K). 

To verify the condition (b) we will use the IMS localization technique adopting the manifold 
version explained in [p^]. 

Using the same functions Jo, Jj,7 as in the previous section, we have on M: 

hi 

We will also identify the functions Jq, Jj.^ with the corresponding multiplication operators. 
Lemma 2.10 (The IMS localization formula). The following operator identity is true: 

(L) H = JoHJo + J2 JhjHJj,j + l[Jo, [Jo,H]] + i ^^[J,, 7, [J,,,, H]] 

in in 
= JqHJq + Jj^^HJj^^ — fj,a^'^\x, dJo{x)) — fj, a^'^\x, dJj^-y{x)), 

where a^^^ is the principal symbol of H, considered as a function on T*M. 
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The proof can be found e.g. in |19]. 

The last two terms in (L) are F-invariant and easily estimated as 0{fi^^^'^), so we have 

(7) H > JqHJo + JjnHJj,i - Cix^'^'^I. 

Since ^Ha > 0, we also have with some constant c > 

(8) JoHJq > fi-^VJl > cfi-'+^^Jl 

if /i € (0, /xo), and is sufficiently small. Note that the coefficient in the right hand side here 
tends to +oo as // | 0. 



Following Lemma 3.4 in |19] we find that 



(9) JjHJj > (1 - Cfi'^)JjKj{fi)Jj - Cfi^^'-^jf . 
But we also have 

Jj^jHJj^^ = Jj^^TjHT^ Jj^^ = T^JjHJjT^ 
> (1 - Cn'')T^JjKj{ij.)JjT-^ - Cfi^''~%jfT-^ 
= (1 - C^,'')T^JjK,{fi)JjT-^ - Cfi^^^'-^Jl^ . 
Let us sum over j, 7 and add D and also JqHJq. Using the inequality 

Kjifi) + Le'^^^ > XI, 
and also the estimates (|7|),(§),(^) above, we obtain then 

H + D> cfi-^+^^J^ + (1 - C/x'=)A J]n " ^/^^""^ 4-1 " Cf^^'^^'I. 

hi 3,1 

Choosing here k = 2/5 we obtain 

(10) H + D>{\-Cfi^'^)I 

with a constant C and with /i G (0, /Uq) for a sufficiently small /io- This proves condition (b) and 
ends the proof of Theorem 2.1. □ 

Remark 2.11. If H is flat near all points Xj, then the estimate (^) can be improved as follows: 
for any R > and e > there exist C > and > such that for all A < and fj, S (0, /xq) 

NriX;H) < N{X + f,^-';K). 

Together with the improved estimate from below this provides the inclusion 

oo 

spec(i/) n (-oo,i2] C ^J(ap-C//l-^ap + Cexp(-C7-V^^^^))• 



The necessary arguments can be found in |19| 
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